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At biome-scale, terrestrial carbon uptake is controlled mainly by weather variability. Observational data 
from a global monitoring network indicate that the sensitivity of terrestrial carbon sequestration to mean 
annual temperature (X) breaks down at a threshold value of 16° C, above which terrestrial C0 2 fluxes are 
controlled by dryness rather than temperature. Here we show that since 1948 warming climate has moved 
the 16° C T latitudinal belt poleward. Land surface area with X > 16 °C and now subject to dryness control 
rather than temperature as the regulator of carbon uptake has increased by 6% and is expected to increase by 
at least another 8% by 2050. Most of the land area subjected to this warming is arid or semiarid with 
ecosystems that are highly vulnerable to drought and land degradation. In areas now dryness-controlled, net 
carbon uptake is —27% lower than in areas in which both temperature and dryness (X < 16 °C) regulate 
plant productivity. This warming-induced extension of dryness-controlled areas may be triggering a 
positive feedback accelerating global warming. Continued increases in land area with X > 16° C has 
implications not only for positive feedback on climate change, but also for ecosystem integrity and land 
cover, particularly for pastoral populations in marginal lands. 



Warming climate is altering climate-control mechanisms of terrestrial carbon sequestration 1-4 . The 
direct observational evidence 5 provided by a global network (FLUXNET) of continuous in situ mea- 
surements of land-atmosphere exchanges of C0 2 , water vapour and energy across biomes and con- 
tinents, indicates that terrestrial C0 2 fluxes are: (1) strongly limited by mean annual temperature (X) of less 
thanl6°C at mid- and high-latitudes; (2) strongly limited by dryness at mid- and low-latitudes; and (3) co-limited 
by both temperature and dryness around the mid-latitudinal belt (45 °N). The sensitivity of terrestrial C0 2 fluxes 
to X breaks down at ~16°C, a threshold value above which no further increase of C0 2 uptake with increasing 
temperature was observed and dryness influence overrules temperature influence. Here, we examine a hypothesis 
that the threshold-latitudinal belt at which Xis 16°C is shifting poleward as the Earth's surface warms and hence 
the areas of dryness-control of terrestrial C0 2 fluxes (X > 16°C) is expanding. We use global land temperature 
data 6 7 (1948-2012) to test this hypothesis and examine the potential consequences of warming-induced exten- 
sion of the dryness-controlled area on climate change. 

We calculated land area where Xis higher than or equal to 16°C for each year during the period from 1948 to 
2012 using mean monthly surface temperature data (0.5° X 0.5° resolution) from the National Centers for 
Environmental Prediction-National Center for Atmospheric Research (NCEP/NCAR) reanalysis dataset. We 
refer to land area with X a 16°C as dryness-controlled areas where terrestrial C0 2 fluxes are limited by dryness 
and not by temperature 5 . A pronounced increase in the dryness-controlled area occurred following a slight drop 
before 1976, mirroring the variation with land warming (Fig. la). About 90% of the variance in the extension of 
the dryness-controlled area was accounted for by land warming (R 2 = 0.90, p < 0.0001, see Fig. lb). 

We assume that net ecosystem- atmosphere exchanges of C0 2 (NEE) in the areas close to the cold side (X < 
16°C) of the shifted boundary are controlled by both temperature and dryness 5 , written as NEE B that is 
predicted by equation (1) (see Methods Summary). For NEE in the area on the warm side (X > 16°C) of the 
shifted boundary, it is written as NEE D , that is determined by dryness alone through equation (2). How will 
increasing X affect NEE in area that shifts from control by both X and dryness (X < 16°C) to control by 
dryness alone (X > 16°C)? We estimated the difference between NEE B and NEE D by applying a NEE model 
(see Methods Summary) that was derived from datasets collected by a worldwide, tower-based, observational 
network 5 to the shifted area. The climate data used in the model were X and dyness averaged over the period 
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Figure 1 | Relationship between dryness-control area (%)and land temperature (°C) (1948-2012): (a) the evolution of dryness-control area (blue line) 
and land average temperature (red line); and (b) correlation between annual dryness-control area and annual land-average temperature (R 2 = 0.90, 
P < 0.0001). The dryness-control area refers to the total area of regions where mean annual temperature was higher than or equal to 16°C and terrestrial 
C0 2 fluxes are controlled by dryness rather than temperature based on the direct observational evidence provided by a global monitoring network 5 . The 
annual dryness-control area and annual land surface temperature during the period from 1948 to 2012 were derived from mean monthly temperature 
data at surface (0.5° X 0.5° resolution) from the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP/ 
NCAR) reanalysis data set 6-7 . The black lines indicate the trends of dryness-control area that was similar to that of land-average temperature (omitted): a 
slight drop between 1948-1975 and then a striking increase during 1976-2012. The striking increase in temperature is a direct result of increased 
greenhouse gases in the atmosphere 31 . The red arrows in (a) indicate El Nino years with oceanic Nino index (ONI) greater than +1.0, while blue arrows in 
(a) indicate La Nina with ONI less than —1.0 (http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml). 70% El Nino 
years were consistent with warmer years, while 80% La Nina years were consistent with cooler years. 



1948-2010 with 0.5° X 0.5° spatial resolution in the shifted area. 
Dryness was calculated from the monthly datasets of net incoming 
short-wave radiation and net long-wave outgoing radiation of the 
NCEP/NCAR reanalysis 8 and the gridded monthly terrestrial pre- 
cipitation datasets 9 . This data-driven estimate indicates that C0 2 
transfer from the atmosphere to the biosphere is reduced by 27% 
in the shifted area where T changed from less than to greater 



thanl6°C. Qualitatively, the model prediction reveals a positive 
feedback mechanism: climate warming extends the dryness-control 
area, which reduces C0 2 transfer from the atmosphere to bio- 
sphere. Thus, because the atmopheric C0 2 concentration will 
increase at a greater rate, the climate will warm at an accelerating 
rate due to the positive feedback. If the global area under dryness- 
control (Fig. la) continues to increase at only the same rate as 



SCIENTIFIC REPORTS | 4:5472 | DOI: 1 0. 1 038/srep05472 



2 




1.0 



(b) Trends 
Wet 



Area below 16 "C 




Figure 2 | The Palmer Drought Severity Index (PDSI). (a) Links between PDSI and ENSO events.The red curve shows PDSI for the area with 
temperature above 16°C; the green curve for the area with temperature below 16"C; the grey curve for the whole area of the land; and the blue curve for the 
shifted area from below 16°C to above 16°C during 1948-2012. The red arrows indicate El Nino years with oceanic Nino index (ONI) greater than +1.0, 
while blue arrows indicate La Nina with ONI less than — 1.0 (http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml). 
(b) The trends of PDSI. The filled circles are five-year moving average of the PDSI data shown in (a). Mean annual land surface temperature during the 
period from 1948 to 2012 was derived from mean monthly temperature data at surface (0.5° X 0.5° resolution) from the National Centers for 
Environmental Prediction-National Center for Atmospheric Research (NCEP/NCAR) reanalysis data set 6-7 . Annual PDSI data were derived from the 
monthly self-calibrated PDSI data (0.5° X 0.5° resolution, spatial range from 60°S to 70°N, http://www.cgd.ucar.edu/cas/catalog/climind/pdsi.html) 11,12 . 
Drought classification by PDSI are: [-0.49, +0.49] -> normal; [-0.5, -0.99] incipient dry spell; [-1.0, -1.99] mild drought; [-2.0, -2.99] -> 
moderate drought; [ — 3.0, —3.99] — * severe drought 10 . The PDSI behaviours to the ENSO events were different between: (1) the area above 16°C (red 
curve in (a)), 90% El Nino years were dryer, while 70% La Nina years were wetter; and (2) the area below 16°C (green curve (a)), 50% El Nino years were 
wetter, while 40% La Nina years were dryer (see Table 1). 



occurred over the preceeding half-century, the warming- induced 
dryness-controlled area will double by 2050. 

With climate warming much of Earth's land has been moderately 
drying since 1976, averged over all land areas, based on annual 
Palmer Drought Severity Index (PDSI) estimates (Fig. 2a) derived 
from the monthly self-calibrated PDSI data 0.5° X 0.5° resolution 



over the spatial range from 60°S to 70°N 10 ". Our analysis finds that 
the drying trend in the shifted area was strongest and in land areas 
where Tis above 16 U C was second strongest (Fig. 2b). The land area 
where T > 16°C encompasses low latitudes in the northern hemi- 
sphere, most of Africa, Middle and South America, Australia, South- 
and Southeast Asia (Fig. 3). In these regions, tower-based FLUXNET 
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Figure 3 | Map of mean annual temperature (1948-2012): below 16°C in light green regions; above 16°C in light red regions, and shift from below 
16°C to above 16°C in purple regions. The map was produced based on the NCEP/NCAR ranalysis data 6-7 (http://www.esrl.noaa.gov/psd/data/gridded/ 
data.ncep.reanalysis.derived.pressure.html), and created using Matlab. The information of vegetation distribution, and precipitation (P) are summarized 
in boxes for the shifted areas (purple color) in each of seven framed regions marked on the map. The vegetation is coded according to the IGBP 
classification: GRA, grassland; CRO, cropland; MF, mixed forest; OSH, open shrubland; WSA, woody savanna; SAV, savanna; EBF, evergreen broad-leaf 
forest; and BAR, Barren or sparsely vegetated. 



observations 12 document that at ground-level these large land areas 
indeed are drying up and this is confirmed by remote sensing data 13 . 
This drying is attributed to increased evaporation and evapotran- 
spiration due to warming. If the trend of drying up over the large land 
area where T > 16°C continues a strong positive feedback on warm- 
ing is suggested because of reduced C0 2 transfer from the atmo- 
sphere to land (expansion of the brown areas in Fig. 3) via NEE 
that is limited substantially by water availability (see Fig. 2b in Ref. 
5), thus inducing additional warming. In contrast, in the land area 
where mean annual temperature is below 16°C (green area in 
Figure 3) a trend toward greater wetness has been observed with 
climate warming (Fig. 2). 

Two large areas between the cold (<16°C, green color in Fig. 3) 
and warm zones (>16°C, brown color in Fig. 3) have different per- 
formances during El Nino/Southern Oscillation (ENSO) events 
(Fig. 2a, Table 1). In this analysis we included El Nino years with 
an oceanic Nino index (ONI) greater than +1.0 during the period 
between 1948 and 2012, and La Nina years with ONI less than — 1.0 
(Table 1). Half of the El Nino years were consistent with the cold 
phases (dips of temperature curve) of the cold part of the land (CPL, 
green area in Fig. 3), while 70% of the La Nina years were consistent 
with the warm phases (peaks of temperature curve) of the CPL 
(Table 1). The CPL warm/cold phases appeared to be opposite of 
the warm/cold phases of the ENSO cycle. However, the global land 
area followed the warm/cold phases of the ENSO cycle very well. The 
CPL wet/dry phases appeared different from that of the global land 
area (Table 1). However, the wet/dry phases of the warm part of the 
land (WPL) were very consistent with that of the ENSO cycle, i. e. 
90% El Nino years were in the dry phases, while 70% La Nina years 



were in the WPL wet phases (Fig. 2a). We could not find a better 
relationship of the WPL warm/cold phases with the ENSO cycle. 
However, if we assume that WPL temperature responses to the 
ENSO cycle lag by a year, 90% of El Nino years coincided with 
the WPL warm phases, while 80% of La Nina years coincided with 
the WPL cold phases. These fascinating coincidences, that became 
obvious after lagging the data by a year, can be understood at least 
theoretically in the following way. In the WPL wet phases, a much 
larger fraction of net radiation is used for evapotranspiration as latent 
heat and hence potential warming is reduced, while in the WPL dry 
phases, comparitively less net radiation is used as latent heat, so the 
temperature is increased. This energy budget adjustment may need 
about a year to reach equilibrium for about half of Earth's land. 
Temperature responses to the ENSO cycle of the shifted area (purple 
color in Fig. 3) were similar to the responses of the total global land 
area because the temperature of the shifted area is close to the land- 
average temperature. However, the pattern with 60% of El Nino years 
being wetter while 30% La Nina years were dryer for the shifted area 
is coincident with the typical precipitation patterns of the ENSO 
cycle reported by NOAA 14 . 

The shifted areas are transitional zones where not only is the 
climate-control mechanism of NEE switched as discussed above, 
but also meteorological conditions are more variable and vegetation 
is highly vunerable to climate changes and weather extremes. 
Dominant vegetations in the shifted regions (purple color in 
Fig. 3) are open shrublands (25%), croplands (22%), grasslands 
(7%), and desert (13%) (see Supplementary Table 1). Except for 
the shifted areas in southeastern China (box 4 in Fig. 3) and south- 
eastern United States (box 1 in Fig. 3), most shifted areas are arid and 
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Table 1 Temperature, 


PDSI, 


and ENSO events from different areas classified in Figure 3 
Properties 1 0 El Nino years* 


10 La Nina years + 


Olobal Land Area 




T ( C) 


70% warmer 


ou /o cooler 






PDSI 


30% dryer 


4\j /o wetter 


Land Area Above 1 6°C 




T (X) 


Not clear but taking 1 year lag 90% warmer 


Not clear but taking 1 year lag 80% cooler 






PDSI 


90% dryer 


70% wetter 


Land Area Below 1 6°C 




T (°C) 


50% cooler 


70% warmer 






PDSI 


50% wetter 


40% dryer 


The shifted area 




T ( C) 


70% warmer 


80% cooler 






PDSI 


60% wetter 


30% dryer 


•10 El Nino years include: 1957-1958, 1965-1966, 1 
+ 1 0 La Nina years include: 1 950-1 95 1 , 1 955-1 956, 


972-1973, 1982-1983, 1986-1987, 1991-1992, 1994-1995, 1997- 
1964-1965, 1970-1971, 1973-1974, 1975-1976, 1988-1989, 1998- 


998, 2002-2003, 2009-2010. 
-1999, 2007-2008, 2010-201 1 . 



semi-arid land with typical vegetation of open shrublands and grass- 
lands. The annual NEE of these ecosystems is quite sensitive to 
climate conditions of low precipitation and high evapo transpiration 
rates 15 - 18 . 

The locations of shifted areas in the northern hemisphere are 
coincident with the descending branch of the Hadley cell (HC) and 
are consequently associated with low precipitation and high evap- 
oration rates 19 . Several lines of evidence indicate that the HC has 
intensified and expanded poleward over the past three decades as a 
consequence of climate warming 20 22 and that the HC expansion in 
the northern hemisphere is stronger than in southern hemisphere 23 . 
The contemporaneous poleward shift of both the HC and the WPL 
(significant since late 1970s) and location of the WPL with respect to 
the HC (nothern desending branch of the HC), strongly suggests that 
the WPL migration poleward is driven by global warming. The dry- 
ing trend of the shifted area (Fig. 2a-b) should be expected to result in 
vegetation cover shift, with decreased biodiversity and desertifica- 
tion. A line of evidence from remote sensing imagery indicates that 
drying is accelerating the degradation of vulnerable shrublands in 
some semiarid Mediterranean area 24-25 . 

Division of the land into the WPL and CPL by threshold value 
(16°C) of annual mean temperature based on 64 years (1948-2012) 
climate data brings new insights into the warming of Earth's surface. 
The two parts of the land behave almost opposite to each other in the 
phases of the ENSO cycle and differ in climate control mechanisms of 
carbon sequestration 5,26 ' 27 . The WPL has expanded poleward signifi- 
cantly (Fig. 1) and has become dryer (Fig. 2) in the past four decades. 
The trend of warming- induced drying of the WPL, by reducing NEE 
thereby reducing withdraw of C0 2 from the atmosphere, contributes 
a positive feedback on global warming. Furthermore, as lands are 
shifted from CPL to WPL becoming more arid and subject to deser- 
tification, they also release soil carbon adding additional C0 2 to the 
atmosphre. It is estimated that 19-29 Pg of carbon were added to the 
atmosphere from vegetation and soil carbon pools globally by deser- 
tification 28 . The frontal boundary (or the shifted area) of the WPL has 
been transformed by global warming into more vunerable regions 
where weather gradients are stronger (Fig. 2), ecosystems are more 
sensitive to even slight increases in water deficit (Fig. 3) 25 , crop yield 
is reduced by extreme heat waves 29 , and vegitated land cover and 
pastroral population are reduced. For instance, in Australia, where 
wide areas are becoming not suitable for sheep breeding due to 
reduced precipitation and increased soil salinity. An expansion of 
the global network 30 monitoring NEE to target the identified shifted 
areas would provide data that could improve our ability both to 
model these regions as they undergo further transitions and to assess 
the likely impacts on climate as a consequence of altered NEE and 
increased soil aridity. The present work raises the following two 
questions: (1) what atmospheric circulation mechanisms support 
the hypothesis of a year time lag between the WPL temperature 
response and the ENSO water phases; and (2) is the synergistic pole- 
ward expansion of the frontal boundary of the WPL with the HC a 



long-term or a short-term behavior and what are the consequences of 
this synergy for global NEE and for the rate of change in atmospheric 
C0 2 ? 

Methods 

Details of calculating land temperature, precipitation, net radiation, and PDSI are 
given in the Supplementary Information. Here we summarize the method used to 
estimate NEE difference induced by the switch of climate control from the CPL to 
WPL. For the case of the shifted area in the CPL (purple in Fig. 3), its NEE B is limited 
by both temperature { T) and dryness (D) and is estimated using a bivariate- nonlinear 
regression model 

NEE B = -3.9855 -0.0272T + 2. 9394D-0. 01 14T 2 -0.69688D 2 , (l) 

where D is defined as RJ (XP), R n is mean annual net radiation MJ m" 2 yr" ', P is mean 
annual precipitation mm yr -1 , and k ( = 2.5 MJ kg -1 ) is the enthalpy of vaporization. 
For another case of the shifted area in the WPL where NEE D is limited by dryness 
along (T > 16°C), we use the regression model of D-limited group in Ref.5, 

NEE D = -1.010LD 3 + 3.1203D 2 + 1.8055D-8.2528. (2) 

All the regression coefficients in equations (l)-(2) are estimated from the published 
eddy-covariance data (see supplementary Table SI in Ref. 5). 
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